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Nonsingular cosmologies are investigated in the framework of f(R) gravity within the first order 
formalism. General conditions for bounces in isotropic and homogeneous cosmology are presented. 
It is shown that only a quadratic curvature correction is needed to predict a bounce in a flat or 
to describe cyclic evolution in a curved dust-filled universe. Formalism for perturbations in these 
models is set up. In the simplest cases, the perturbations diverge at the turnover. Conditions to 
obtain smooth evolution are derived. 
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I. INTRODUCTION 

Singularities occuring in General Relativity (GR) may 
be avoided in more fundamental frameworks where GR 
predictions are recovered as the low energy limit. Such 
frameworks could be provided by String or M-theory, or 
loop quantum gravity. In particular, the evolution of the 
Big Bang cosmology might be extended to a preceding 
contracting phase which, due to new physics relevant at 
high curvature or energy scales, turns into the expanding 
phase our universe is experiencing now [l|, 0| , thus avoid- 
ing the Big Bang singularity (which in the inflationary 
picture [3| is manifest rather as geodesic incompleteness 
[J] than divergence of curvature invariants). These sce- 
narios are called bouncing cosmologies @. 

Bouncing cosmologies can solve the horizon problem, 
but to replace inflation they should, among other things, 
also predict a viable, nearly scale invariant spectrum of 
perturbations. This tricky issue can be circumvented if a 
curvaton field is responsible for the generation of fluctu- 
ations 6]. Otherwise, matching conditions are often re- 
quired to track the evolution of the perturbations across 
the bounce (tJ. It has been noted that the curvature 
perturbation may become singular at the bounce, while 
the gravitational potential, whose growing mode usually 
persists in the post-bounce era, may have regular behav- 
ior Q . A general solution for the perturbations supports 
these conclusions It is also well known that the fea- 
tures of the spectrum can depend sensitively upon the 
details of the dynamics of the bounce and the physics 
behind it. 

Hence it is useful to consider explicit examples which 
allow one to scrutinize the possible behaviors of fluctua- 
tions at the bounce. However, it is necessary to violate 
energy conditions (EC). At least the strong EC must be 
broken to change the sign of the expansion rate, and the 
null EC cannot be respected if there is no curvature. This 
rather generically introduces pathologies that, though 
one may interpret them as only a shortcoming of the ef- 
fective theory, hinder from reaching definite conclusions 



of the evolution of the spectra 1 1 Ol — il 2f | . An example is the 
perturbation divergence in pre-Big Bang cosmology [l3| 
which can be shown to be an indication of an appearance 
of a ghost [14|. To avoid EC violating matter fields, one 
can contemplate modifications of gravity that introduce 
no new degrees of freedom. This can be achieved with 
an action involving an infinite series of d'Alembertians 
acting on the curvature invariants in such a way that the 
propagator has no poles; these string-inspired nonlocal 
modelfo have been shown to be ghost-free and asymp- 
totically free at their Newtonian limit [IH [26| . 

In the present paper we consider the simpler case of 
second order f(R) gravity, which has similar desider- 
able features. The fourth order metric f(R) models cor- 
respond to scalar tensor theories of the form C^th = 
<f>R + V ((/)). By erasing the kinetic term implicit in the 
nonminimal coupling, one obtains the second order the- 
ory C 2n d = 4>R- j^{d(t)) 2 + V{4>). Though this can be 
problematic in view of the well-posedness of the Cauchy 
problem [27j . these simple models may avo id g eneric in- 
stabilities present in higher order theories [28j and thus 
provide an effective description of low energy effects of 
quantum gravity. In particular, loop quantum gravity is 
expected to modify the cosmological dynamics at high 
curvatures without introducing new degrees of freedom. 
To obtain the quadratic density correction appearing in 
the particular scalar loop quantum cosmology, one may 
need to consider an infinite number of terms in the po- 
tential V(<p), which may be interpreted to reflect the 
nonlocal nature of the underlying theory as discussed in 
Ref . . Note also the recent extension of the framework 

Low curvature corrections in these so called Palatini- 
f(R) theories have also been considered as alternatives 
to dark energy [3ll l32j , but though they may produce vi- 
able background expansion, they generically fail to pro- 
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1 A biscalar-tensor model [la . motivated by the nonlocal cos- 
mology based on inverse dAlembertian gravity [l7l. 1 1 SH may also 
accommodate bounces [19j |. though their viability remains to be 
shown 20 1 . Bounces in modified gravity were considered also in 
e.g. [IJEil. 
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duce the observed matter power spectrum, at least for 
pressureless dust cosmology |33l - [35| . Problems may ap- 
pear also at microscopic level, as discussions of electron- 
electron scatterin g an d Hydrogen atoms, seem to imply 
[3ol \st\ , see also |38l , |39| . This may be due to need to 
reconsider the averaging problem in these models pol |. 
or the coupling of gravity to matter taking torsion and 
nonmetricity into account [4l|-|43j. One may adopt the 
approach of considering the formalism as an effective 
macroscopic description, and then new phenomenology 
can emerge from the potentially viable high curvature 
corrections. Studies of spherically symmetric systems 
show that the classic Solar system tests are passed by 
these models while the high curvature effects have 
interesting predictions for white dwarfs and neutron stars 
[45j . Bounces have been suggested too |46| . 

We briefly review the first order formalism approach to 
nonlinear curvature gravity in section [TT] where we also 
write and solve the cosmological background equations 
taking into account spatial curvature. We derive the con- 
ditions for bounces to occur and confirm them numeri- 
cally. In section IIIII we present the equations governing 
the evolution of perturbations in convenient forms, based 
on derivations in Refs.j47l.l48l]. Corrections to and gener- 
alizations of previous literature are pointed for both the 
background and the fluctuation equations. In section ITVl 
we discuss these results and their implications. 



II. BOUNCING BACKGROUNDS IN 
PALATINI-F(R) GRAVITY 

After writing the general equations for the generalized 
gravity model, we derive the bouncing conditions and 
analyze them both analytically and numerically. 



A. Palatini approach to generalized gravity 



Consider gravity theories represented by the action 



S = d n x 



(1) 



respect to f 2 , obtains 



j^F] =0. 



(4) 



implying that this connection is compatible with the con- 
formal metric 



g M „ = F 2 ^ 



g^ ■ 



(5) 



This connection governs how the tensor R^ v appearing in 
the action settles itself, but it turns out that the metric 
connection determines the geodesies that freely falling 
particles follow, since the energy momentum 



_ 2 6(J=g£ m ) 
is conserved according to this connection, 



(6) 



(7) 



whereas in general V M T^ ^ 0. Therefore we have a met- 
ric theory of gravity [49| in the sense of Ref . [50l] . The 
trace of the field equations allows us to solve R as an 
algebraic function of the matter trace T = g^T^. This 
central relation reads 



FR — 2f = T. 



(8) 



Prom now on we set the spacetime dimension to n = 4 
and use units 8ttG = c = 1. Written in the form of GR 
plus correction terms, the field equations read: 



(l-F)RZ(g)-—(V»F)(V v F) 



(/-fl) + (l--)DF- 



2F 



(dFf 



Since the corrections can be expressed as functions of 
the matter trace, one can view Eq.© as GR with gen- 
eralised coupling to matter: only the way that "matter 
tells spacetime how to curve" is modified. So, the whole 
RHS may be regarded as an effective matter energy- 
momentum tensor. In vacuum it reduces to a cosmo- 
logical constant [•[> lj . This is also the case in the presence 
of conformal matter, i.e. if T = 0. 



Here (/>,... are some matter fields. In the Palatini ap- 
proach one lets the torsionless connection r^j vary inde- 
pendently of the metric. The Ricci tensor is constructed 
solely from this connection, 



Ruu 



-pa -p\ 



f ° A f x av . (2) 



The field equations which follow from extremization of 
the action Eq.([T]) with respect to metric variations, can 
be written as 



FR?, 



(3) 



B. Background cosmology 

In the spatially flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) universe with the line element 



ds 2 = -dt z + a 2 (t) 



dr 2 



1 - Kr 2 



r 2 dY? 



(9) 



and a perfect fluid source with a constant equation of 
state w = p/ p, the Friedmann equation can be written 
as 



where we have defined F = df/dR. In GR, (R — 
2A)/8"7rG, so F = l/8wG. By varying the action with 



6F ( H+ 2Fh) =P + 3P+f-S F §- (10) 
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In the Appendix [X] we write the general Friedmann equa- 
tion when w = is not assumed. The Hubble parameter 
can be expressed fully in terms of the curvature scalar, 
and in our case then rewritten as 



H 2 = 



1 3(1 + w)f - (1 + 3w)FR- 6Ff 



K 



(l-3w)F 



~ 2K L+W ) F(BF'-F) J 



(11) 

We have expressed the Hubble rate as a function of R, 
which we in turn may solve from the trace equation (|8]l. 
If the scale factor is monotonic, one may find its evolu- 
tion algebraically once the given matter content as out- 
lined in Ref. [52[. In our case however it is preferable to 
solve the numerical system by the integrating differential 
equations, which are shown in the appendix. 



C. Bouncing solutions 

In general, a necessary condition for a bounce to occur 
is obtained from (fTT|) and j8]) as 



F = 0. 



or 



K 



F(l2—-R)=3p(l + w). 



(12) 
(13) 



Now, by solving the trace equation ([5J, FR is given by 
an inverse function of the density, and its form is very 
dependent of f(R). The quadratic term, which may be 
considered the leading correction^ to GR, can already 
lead to bouncing cosmology. Let us thus consider the 
case 



f(R) = R + aR 2 , 



(14) 



to study explicitly the background behavior. This model 
results in a symmetric bounce (the pre-Big Bang is the 
time reversal of the post-Big Bang evolution). The Fried- 
mann equation is now 



2a 3 (a 3 - aRo)'' 



■mi- : ., ., >/ :i /A, - %Ka{a 3 + 2,.//,,; + a Rl] 

(15) 

where Rq > is a constant which is equal to the matter 
density at a = 1. The derivatives of the Hubble rate are 

~^U). The bounce 



(16) 



written in the appendix as (|A3I) and 
condition (|1 21 13[) now becomes 



a 3 + 2aR 
(2a 3 + aR )R 



, or 
6Ka(a 3 



2aR Q ) 



In the flat case, the second condition is simply a 3 = 
— aRo/2, which can be satisfied given a negative a < 0. 



2 In addition, one notes that this is the only case when the trace 
equation JS) is linear in the sources. One may then contemplate if 
more general functions f(R) would be, after averaging, effectively 
described by the quadratic model. 



However, the first condition will be saturated earlier 
when the scale factor has contracted to a 3 — — 2ai?o. 
Then, at the bounce we have F —> 0. 

Let us consider whether one may avoid F = at the 
turnover in curved models with K =/= 0. Since the nor- 
malization of the scale factor is arbitrary, let us assume 
here that the second bounce condition is fulfilled before 
the first one at a = 1. This implies that 



2aR > -1 



and that 



R = 6K 



1 + 2aR 

2 + aR 



(17) 



(18) 



Solving a from the second constraint and plugging into 
the first condition gives 



4(3K-Rp) 
Ro - 12K 



> -1 



(19) 



We should assume Rq > 12K, since in any realistic uni- 
verse the curvature is subdominant to the matter density 
at early times by many orders of magnitude. Then the 
constraint reduces to Rq < 0, in contradiction to our as- 
sumptions. If we however allow a positive curvature to 
dominate over matter density at the bounce, K > R /12, 
we can realize bounces where F stays finite at the turning 
point. This could have occurred if there was significant 
amount of inflation after the bounce which diluted away 
both the curvature and matter density. However, the de- 
tails of such a case are not of interest to us, as there also 
the possible signatures from bounce were most probably 
erased. 

As we will see in the following, the F — bounce re- 
sults in divergence of perturbations. In general, assuming 
negligible curvature but allowing general f(R) and w, the 
condition for the second type of bounce (TIB"]) assumes the 
very simple form 



F>0, f(R) = -p(l + 3w) 



(20) 



In the simplest models considered here in detail, this is 
a necessary condition for perturbations to stay regular 
at the bounce. In section HVl we briefly discuss possibly 
viable generalizations of the models. 



III. PERTURBATIONS 

In the following we first specify our perturbation sys- 
tem, then consider the evolution of the perturbations in 
pressureless matter as a specific example, and its descrip- 
tion in terms of the canonical variable. Cosmological 
perturbation theory is presented in the reviews (H3, |54|. 
and applied to this class of generalized gravity theories 
inRefs. 1,1. 
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FIG. 1. Solid (black) line: cyclic evolution of the scale factor 
in the quadratic model with K > 0. Dashed (blue) line: the 
same model with doubly as much matter. Dash-dotted (red) 
line: the same model with tripled matter density. 



A. The perturbation system 

The line-element in the perturbed Fricdmann- 
Lemaitre-Robertson- Walker (FLRW) spacetime can be 
written as 



ds 2 



- (1 + 20) dt 2 



(3) 
Hj 



hi 



b l dx l dt 



(21) 



We characterize the scalar perturbations in the longitu- 
dinal Newtonian gauge by the variables gravitational po- 
tentials 4> and ip. Vector perturbations introduce two 
more degrees of freedom, encoded here into the diver- 
genceless 3-vector field hi. Gravitational waves are de- 
scribed by the two free components of the symmetric, 
transverse and traceless 3-tensor hij . The comoving spa- 

(3) 

tial background metric gr>. reduces to in a flat uni- 
verse. The vertical bar indicates a covariant derivative 

(3) 

based on the Levi-Civita connection of q). ■ . This metric 
is used to lower and raise spatial indices i,j,k... of the 
perturbation variables. The components of the energy- 
momentum tensor for a general fluid is imperfect fluid 
are 



rpO 

1 o 



-(p + Sp), 



T? = -(p + p) (v ti 



(22) 
(23) 

T i j ={p + 6p)5)+W j . (24) 

Here p and p are energy density and pressure, and v, 
are the scalar and vector velocity perturbations, respec- 
tively. Background quantities are denoted with an over- 
bar, which we will usually omit when unnecessary. The 
isotropy of the background does not allow anisotropic 
stress except as a perturbation. This we decompose into 
the scalar, vector and tensor contributions as 



Hi 



3 



n 



n 



(*) 



(25) 



where A stands for the three-space Laplacian based on 



f 31 i ■ 

the Levi-Civita connection of g\j . The vector 11^ is 

divergence-free and the tensor IL[*^ is symmetric, trans- 
verse, and traceless. This completes our specification of 
the perturbation system. 



t(v) 



B. Fluid quantities 

Next we will discuss the evolution of the system in 
terms of fluid variables. First we consider the density 
perturbation in the comoving gauge (i.e. CDM rest frame 
in the present case) and then the velocity perturbation 
in the uniform-density gauge (i.e. the frame where CDM 
is smoothly distributed). The former quantity becomes 
ill-defined at the bounce, the latter behaves somewhat 
better. 

It is convenient to introduce the comoving density per- 
turbation A which is given by longitudinal gauge quan- 
tities as follows: 



A = 6 



3H(l + w)^ 



(26) 



The evolution equation for the perturbations has been 
derived in the general case and is of the form 



A = D X HA 



k 1 

D 2 H 2 + D k —)A + Pxffn + P 2 H 2 U , 



(27) 

where the dimensionless coefficients are given in the ap- 
pendix of Ref. [HI . In the case of pressureless dust the 
evolution equation simplifies to 



H 



H 



k 2 



A + (2H + T)A=[- + 2H+-T-^ci ff )A 



H 



where we have defined the auxiliary quantity 
2 



(28) 




and the effective sound speed squared 



F 



L eff 



Z(F + 2FH) 



(30) 



If both T and c 2 ^ vanish GR evolution is recovered, 
so these variables represent the modified gravity effects. 
In the specific example model discussed in section III C[ 
both of these terms apparently diverge at the bounce. 
In particular, we have divisions by H, F + 2FH and F, 
where the first term vanishes always, the second term 
vanishes at least for dust, and the last term vanishes at 
least for the flat dust bounce model. In particular, for 
the quadratic model 



4 



a 3 — aRo 



(31) 
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T = 



ISHaRoa 3 



(a 3 - <xR )(a 3 + 2aR a ) ' 



(32) 



and we see that when a = —aRo, the sound speed in 
fact is regular but the factor T is not. The H and H/H 
are given in the appendix as (|A3[) and (|A4I) . Though the 
comoving gauge is where we want our observable density 
in, this coordinate system can become ill-defined at the 
bounce. However, a possibility remains that this is not a 
physical problem. 

The perturbations can be carried across the bounce in 
another variable than the comoving gauge density per- 
turbation. Another convenient quantity to consider is 
vg, the velocity perturbation of matter evaluated in the 
uniform-density gauge. From Eq. (l26[) . we have 



3aH 



-A. 



(33) 



and readily obtain from (|27|) the evolution equation 
which coincides with the Eq.(46) in Ref . [47l| . 

v s + [AH + 2— + f\ v s 



+ I 3(if + H 2 ) + HJ r +-^c 2 ff J v$ =0,. (34) 

The divergent 1 / H terms in the second expression cancel, 
H 



AH + 2'—+F 
ti 



F + 2FH 



F 2 

F + 2(FH)' +AFH 2 

F 



.(35) 



However, the F/F can divergent, when F — > in the 
flat bounce model. This can be avoided with curvature 
K > p/12, or in the bounce of the type (|2H)) but even 
then typically F —¥ at the turnover. Consequently, 
the factor T remains apparently divergent due to T ~ 
Thus, by considering (|3"31 instead of (|2"rJj) the 



F+2FH 

apparent divergencies of the coefficients in the evolution 
equations can be made less severe (from ~ 1 /H 2 to ~ H ) 
but they seem to persist. 



Canonical variable 



The canonical variable v obeys the equation of motion 
i)+Hu+ (--H 2 + Ih-*+C„+ k —^c 2 ef ^\ v = , 



where 



4 2 a 2 



c » = -&( 3H+ f' 



(36) 



(37) 



The coefficients remain apparently divergent for the flat 
dust bounce model of section III C[ because the potential 



becomes infinite due to the ~ F 2 / F 2 term. However, for 
the curvature-dominated model or bounce of the type 
([20)) both C„ and the sound speed remain regular at the 
bounce, since H,F —> 0. Thus we have found the explicit 
conditions for the perturbations to be carried smoothly 
across the bounce. The relation between the Mukhanov- 
Sasaki variable v and the comoving density perturbation 
is given by 



u(t) = exp 



[H(t') + T(t')} dt' \ 8(t) 



(38) 



Hence, for a model with F not dipping to zero, one may 
use the regular equation (f36]) to solve the evolution, and 
the transformation ([3"B"| to obtain the results in terms of 
the observables. 



IV. CONCLUSIONS 

Within the Palatini framework one may consider ex- 
tensions of GR without introducing new degrees of free- 
dom. Hence they may serve as useful toy models describ- 
ing more completely the cosmological evolution, with mo- 
tivations from e.g. loop quantum cosmology. We showed 
that there are nonsingular bouncing backgrounds in sim- 
ple examples of such quantum corrected gravity models, 
and set up the formalism for the perturbations in these 
models in order to monitor the evolution of their spectra 
across a bounce. 

The models of the type (|12"j) . characterized by F — >• 
at the bounce, were found to feature singular behavior 
of perturbations in a flat, dust-filled universe. This may 
be cured in the curvature-dominated case, reflecting the 
fact that in the K = case the effective matter sources 
necessarily violate the null EC, whereas this is not the 
case if curvature is present. Since F gives the sign of 
the graviton action, pathologies were to be expected at 
F — » 0. At this point the conformal relation between 
Einstein and Jordan frames (JSJ becomes ill-defined. Let 
us note though that as the perturbations explode, their 
backreaction will render the perturbative system invalid, 
and as these nonlinear effects are very difficult to tackle 
in practice, we cannot say if there is a true singularity or 
not and whether the bounce occurs or not. In this light 
the problem is rather unpredictivity. 

Obviously, it would be interesting to explore possible 
ways to obtain bounces of the type (|20)) . Such might be 
constructed by considering just more general functions 
f(R) than the one involving solely a monomial correction. 

Another way to obtain smooth evolution could be to in- 
clude more general sources than completely pressureless 
fluids. Apart from allowing w ^ 0, one may also con- 
sider stabilizing the system with entropic or anisotropic 
stresses. Indeed, this has previously proved successful in 
eliminating instabilities in matter perturbations in some 
dark energy models [48|,[55[ (however, those models based 
on infrared gravity modifications may be otherwise prob- 
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lematical as discussed in the introduction). More realis- 
tically, also radiation would be included as a source and 
this changes the dynamics and possibly the conclusions. 
Our results can be directly applied to such more general 
models with possibly regular evolution. This is left for 
future studies. 

As a concluding remark we note it is possible that a 
complete evolution of the background and structures of 
the universe is not amenable to classical description by 
second order differential equations, and one may have to 
take into account in a more nontrivial way the presently 
unknown physics at the very high curvature scales in or- 
der to provide a fully consistent coarse-grained picture 
of the cosmology that emerges. Meanwhile, the quest for 
the effective field equations for gravitational interactions, 
at both high and low curvature regimes, is ongoing. 



Appendix A: Friedmann equations 
1. General case 



9wF' + (F- F'R)^J3F (F(R - 12§) + 3p(l + to)) 
3{2F 2 -2FF'R-3F'pa w ) ' 

where a w = (1 + w)(l - 3to) and F' = dF/dR. For a 
constant equation of state w = 0, this reduces to 



3Jj2 = £[Mi + ^)-^(i2f -R)] 



F 



3 F> p(l+w)(l-3w) 
2 RF' — F 



(A2) 



In the limit K = 0, this agrees with the formulas in the 
references [1^, H3, [52[ (and in several later references). 
However, this does not reduce to the Hubble law con- 
sidered in Ref.[46|, and consequently our solutions and 
results in subsection III CI are somewhat different from 
theirs. For example, their Eq.(20) for the specific case of 
the quadratic model is quite different from our Eq. (|15|) 
also when K = 0. Note that in Ref.[56| the authors have 
corrected a mistake regarding nonzero spatial curvature 
in Ref . (4|| , and there the considerations are in accordance 
with ours here. 



Without assuming a constant equation of state or van- 
ishing curvature, the Hubble parameter may be written 
as 



H 



(Al) 



2. Quadratic model 

The first two derivatives of the Hubble rate in the 
quadratic model (fl4| are 



H 



a 7 R Q (a 2 - 2AaK) - 2a w K + 6a 4 aR 2 (a 2 - 6aK) + aa 2 R% (3a 2 + 8aK) - a 3 i? 4 

2 (a 4 - aaRo f 



(A3) 



H 



.10 



R (3a 2 - IMaK) - 4a 13 K 



H 2 (a 4 - aaRof 

+ 6a 7 aRl (7a 2 - MaK) + a 4 a 2 i?jj (45a 2 + 52aK) - 4aa 3 i?^ (3a 2 + AaK) + 3a 4 i?| 



(A4) 



Ea. (IA3l) is a suitable form for numerical integration. We 
checked that the solution gives Eq. (fT5|) and that its nu- 
merical derivative reproduces Eq. (|A4p . The latter al- 
gebraic expression is needed specifically at the bounce 
where H, H — > 0. 



ACKNOWLEDGMENTS 

The author is grateful to T. Biswas and A. Mazumdar 
for enlightening discussions and to G. Olmo, P. Peter 
and T. Sotiriou for useful comments on the manuscript. 
This work was supported by the FOM and the Finnish 
Academy. 



[1] G. Veneziano, |Phys. LetEl B265 287 (1991). 

[2] M. Gasperini and G. Veneziano, Astropa rt. Phys.] l, 317 

(1993), [arXiv:hep-th/9211021| 

[3] A. A. Starobinsky, |Phys. LettT| B91, 99 (1980^ 

[4] A. Borde, A. H. G uth, and A. Vilenkin,|Phys. Rev. Lett.| 

90, 151301 (2003), [a^Crv:gr-qc/01106T2r ~ 



[5] M. No vello and S. E. P . Bergliaffa, [Phys. Rept.| 463, 127 
(2008) , |arXiv:0802.1634 [astro-ph]| 

[6] K. Enq vist and M. S. Sloth |Nucl. PhysT] B626, 395 
(2002) , |arXiv:hep-ph/0109214| 

[7] J. Martin and P. Peter , |Phys. RevT] D69, 107301 (2004), 
|arXiv:hep-th/0403173"l 



7 



[15] S. Nojiri and S. D. Odintsov, 


Phys. Lett. B659, 821 [39] 


(2008), arXiv:0708.0924 [hep-th] 





R. Du rrer and F. Vernizzi, [ Phys. Rev.| D66, 083503 
(2002), |arXiv:hep-ph/0203275[ ~ 

V. Bozza, JCA P 0602, 009 (2006), 
|arXiv:hep-th/0512066| 

Tsujikawa, R. Brandenberger, and F. Finelli, 
|Phys. Rev.| D66, 083513 (2002), |arX"iv:hep-th/0207228| 
J. Martin, P. Peter, N. Pinto Neto, and D. J. Schwarz, 
|Phys Rev. D65, 123513 (2002), [arX"iv:hep-th/0112128| 
R. Kallosh, J. U. Kang, A. D. Linde, and V. Mukhanov, 



JCAP 0804, 018 (2008), |arXTv:0712.2040 [hep-th]| 
S. Kawai, M.-a. Sakagami, and J. Soda, |Phys. LetE] 
B437, 284 (1998), |arXiv:gr-qc /9802033 
T. Koi visto and D. F. Mota, |Phys. Rev] D75, 023518 



T. Koivisto, |Pbys. Rev.| D77, 123513 (2008), 
|arXiv:0803.3399 [gr-qc] | 

C. Wetterich. IGen. Rel. Grav.l 30, 159 (1998), 
|arXiv:gr-qc/9704052| 

S. Dese r and R. P. Woodard,|Phys. Rev. Lett.| 99, 111301 
(2007), |arXiv:0706.2151 [astro-phJ ~ 

H. Farajollahi and F. Milani(2010), 
|arXiv:1004.3512 [gr-qc]] 

T. S. Koivisto, |Phys. Rev.| D78, 123505 (2008), 
|arXiv:0807 3778 [gr-qc]] ~ 

I. Y. Aref'eva, L. V. Joukovskaya, and S. Y. Vernov, 

IJHEPI 07. 087 (2007), |arX~iv:hep-th/0701184| 

T. Saidov and A. Zhuk(2010), |arXiv:1002.4138 [hep-th]] 
Y.-F. Cai and E. N. Sa ridakis, JCAP 0910, 020 (2009), 
|arXiv:0906.1789 [hep-th"]j 

L. R. Abramo, I. Yasuda, and P. Peter, P hys. Rev.| D81, 
023511 (2010), [arXiv:0910.3422 [hep-th] | 
T. Biswas, A. Mazumdar, and W . Siegel, JCAP 0603, 
009 (2006), |arXTv:hep-th/0508i"94"l 

T. Biswas, R. Brandenberger, A. Mazum- 
dar, and W. Siege l, IJCAPI 0712, 011 (2007), 
|arXiv:hep-th/0610274] 

N. Lanahan-Tremblay and V. Faraoni, 
[ Class. Quant. Grav.| 24, 5667 (2007), 

|arXiv:0709.4414 [gr-qc]] 

T. S. Koivisto(2009), |arXiv:0910.4097 [gr-qc]| 

G. J. Olmo and P. S ingh, IJCAPI 0901. 030 (2009), 

|arXiv:0806.2783 [gr-qc]] 

L. Fatibene, M. Ferraris, and M. Francaviglia(2010), 
|arXiv: 1003. 1619 [gr-qc] | 

D. N. Vollick, |Phys. Rev.| 
|arXiv:astro-ph/030 6630 
G~ Allemandi, AT Borowiec, 
S. D. Odintsov, |Phys. Rev.| 
|arXiv:gr-qc/0504057| 



D68, 063510 (2003), 

M. Francaviglia, and 
D72, 063505 (2005), 



[33 
[34 
[35 

[36 
[37 
[38 



[40 

[41 

[42 

[43 

[44 

[45 
[46 

[47 

[48 

[49 

[50 

[51 

[52 

[53 
[54 
[55; 
[56 



T. Koivisto, |Phys. Rev.| D73, 083517 (2006), 
arXiv:astro-ph/0602031 



K 
Class 



B. Li, K. C. Chan, and M. C. Chu, |Phys. Rev . D76, 
024002 (2007), [arXTv:astro-ph/0610794| 

Uddin, J. E. Lidsey, and R. Tavakol, 
Quant. GravT] 24, 3951 (2007), 

niXiv: 0705T0232 [gr-qc]| 
Flanagan, 



E. E. 



|Phys Rev. Lett.| 92, 071101 (2004), 
|arXiv: astro-ph /03081 1 1 1 

G. J. Olmo, |Phys RevJ, D77, 084021 (2008), 
|arXiv:0802.4038 [gr-qc] | ~ 
A. Iglesias, W. Kaloper 
M. Park, |Phys. Rev.l D76, 
|arXiv:0708.1163 [astro-ph]j 

Sotiriou 
25, 



A. 



Padilla, and 
104001 (2007), 



and 
062001 



J. C. 

(2008), 



E. Barauss e, T. P. 
Miller, Class. Qu ant. Grav.| 
|arXiv:gr-qc/0703132| 
B. Li, D. F. Mota, and D. J. Shaw, IPhys. Rev.| D78, 
064018 (2008), [arXTv:0805.3428 [gr-qc] [ 
T. P. Sotiriou aiid S. Liber ati, |Annals Phys.| 322, 935 
(2007) , |arXiv:gr-qc/0604006 
S. 



Capozziello 



and 



Vignolo(2009), 



|arXiv:0910.5230 [gr-qc] I 

T. P. Sotiriou, |Class. Quant. Grav.| 26, 152001 (2009), 
|arXiv:0904.2774 [gr-qcTT 

K. Kainulainen, V. Reijonen, and D. Sunhede, 
|Phys. Rev.| D76, 43503 (2007), |arXiv:gr-qc/0611132| 
Reijonen(2009), [arXiv:0912.0 825 
Barragan, G. J. Olmo, and 

D80, 024016 (2009), [arXiv:0907 0318 [gr-qc]| 
and H. Kurki-Suonio, Class. Quant. Grav. 



V. 
C. 

Phys. Rev. 
T. Koivisto 
23, 2355 (2006), 
T. Koivisto, 



H 



gr-qc] 

Sanchis-Alepuz, 



|arXiv:astro -ph/0509422 
IPhys. Rev I D76, 043527 (2007), 

|arXiv:0706.0974 [astro-ph] | 

T. Koivisto, |Class. Quant. GravT| 23, 4289 (2006), 
|arXiv:gr-qc/0505128| 

C. M. Will, Liv ing Rev. Rel. 9, 3 (2005), 
arXiv:gr-qc/0510072 



M. Ferraris, 



Class. Quant. 



M. 
Grav. 



Francaviglia, 
11, 



arXiv:gr-qc/9303007 
M. Amarzguioui, O. Elgaroy, 



T. Multamaki, 


Astron. Astrophys. 


arXiv:astro-ph/0510519 





and 
1505 

D. F. 
454, 



I. Volovich, 
(1994), 

Mota, and 
707 (2006), 



H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 
78, 1 (1984). 

V. F. Mukhanov, H. A. Feldman, and R. H. Branden- 
215, 203 ( 1992). 

F. Mota, |Phys. Rev.| D73, 083502 



berger, |Phys. Rept. 
T. Koivisto and D. 
(2006), ||arXiv:astro-ph/0512135| 
C. Barragan, G. J. Olmo, and H. 



arXiv: 1002.3919 [gr-qc 



Sanchis-Alepuz(2010), 



